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ABSTRACT 

A simple algorithm for the in-flight determination of the magnetic bias of a spacecraft is presented. 
The algorithm, developed for use during the Hubble Space Telescope mission, determines this bias 
independently of any attitude estimates and requires no spacecraft sensor data other than that from 
the spacecraft magnetometer(s). Estimates of the algorithm's accuracy and results from a number 
of numerical studies on the use of this algorithm are also presented. 

I. INTRODUCTION 

I 

The determination of the magnetic bias of a spacecraft before launch can be quite difficult, 
particularly if the spacecraft is fairly large. The lack of such a bias determination can severely 
restrict the usefulness of onboard magnetometers in measuring the ambient magnetic field and thus 

i 

also in determining spacecraft attitude. Attitudes based on magnetometer data are often the first to 
be determined during a spacecraft’s mission. Having an algorithm that can determine the spacecraft 
magnetic bias at the location of each magnetometer without requiring knowledge of the spacecraft 
attitude is therefore desirable. We present here a simple algorithm that fulfills this need. The 
algorithm, developed specifically for the Hubble Space Telescope (HST) mission, requires no 
spacecraft sensor data other than that from the spacecraft magnetometer(s). The only other required 
input data are magnetic field values from a data base geomagnetic field model. Section II of this 
article presents the algorithm; Section III, an analytical error model; Section IV, numerical examples 
of the use of the algorithm; and Section V, our conclusions. 
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II. MAGNETIC BIAS DETERMINATION ALGORITHM 


The spacecraft magnetic bias, D, at the location of an onboard magnetometer is related to the 
measured magnetic field, M, and the local geomagnetic field, B, as follows: 

M - D = A B , (1) 

where M and D are defined in spacecraft coordinates and B in geocentric coordinates, and A is the 
(unknown) spacecraft attitude transformation matrix. The bias D is assumed to be time indepen- 
dent in the spacecraft reference frame. The value of B is to be obtained from a geomagnetic field 
model. Obtaining B from the model requires both a spacecraft ephemeris and the magnetometer 
measurement times. The final accuracy of the determination of D will depend on the accuracies of 
the magnetometer calibration and the geomagnetic field model (including differential errors intro- 
duced by spacecraft ephemeris errors and timing errors), as well as noise in the magnetometer 
readings. Taking the scalar product of each side of equation (1) with its transpose (noting that 
A T = A' 1 ) results in 

M 2 - 2 M T D + D 2 = B 2 (2) 

Equation (2) is independent of spacecraft attitude because B 2 is invariant under rotation. We have, 
however, reduced a three-component vector equation to a scalar equation. Because D has three 
components, measurements must be taken in such a way as to provide at least three linearly 
independent versions of equation (2). This can be done either by taking measurements at a number 
of positions in the orbit or by changing the spacecraft attitude. With the latter approach, it is still 
not necessary to know what the various attitudes are — only that they are distinct from one another. 

Equation (2) may be used as the basis for constructing a least-squares loss function for use in 
solving for D. The loss function is constructed using an error variable, ej, defined by 

ej = (Mj 2 - 2Mj T D + D 2 ) - Bj 2 , (3) 
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where the subscript i indicates measurement number. The corresponding least-squares loss 
function, J, is given by 


J = 1/2 £ Wj ( e i - < ej > ) 2 (4a) 

= l/2Zw i [(M i 2 - 2 Mj T D + D 2 ) - Bj 2 - < ej > ] 2 , (4b) 

where w ; is a weight, inversely proportional to the variance of e i; associated with the ith measure- 
ment, and < ej > is the expectation value of e^ The implicit assumption made in selecting equation 

(4a) as the loss function is that each ej is an independent random variable with Gaussian (or nearly 
Gaussian) distribution. If this assumption is true, then the value of D that minimizes J is the 
optimal estimate for the true spacecraft bias. Specification of the values for w ; and < > is not 

necessarily a trivial matter. The simplest approach, of course, is to assume w ; = 1 and < e ; > = 0 
for all i. We present in Section in a simple analytic model from which more exact expressions for 
W; and < ej > can be derived. It is worth bearing in mind, of course, that because of the uncertain- 
ties involved in modeling errors (both in their magnitude and the form of their distribution), the 
concept of an "optimal'' solution for D becomes a bit fuzzy. All solutions for D within a small 
error window about the value that minimizes I will be essentially equally valid. 

We now construct an algorithm to determine the value of D that minimizes J. Equation (3) is 
quadratic in the unknowns, i.e., in the components of D. This is awkward because it prevents 
simply applying the standard tools of linear algebra to construct a least-squares solution for the 
problem. We note, however, that if the quantity D 2 on the right-hand side of equation (3) were a 
constant, then standard linear techniques would be applicable. We therefore construct a modified 
version, J c , of the loss function: 


J c = 1/2 £ Wj [ ( Mj 2 - 2 MT D + D c 2 ) - B? - < ei > ] 2 , 


(5) 
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where the term D c 2 is to be considered an adjustable parameter rather than a function of the 
components of D. Standard linear least-squares techniques (e.g., see Reference 1) may be 
used to find the value of D that minimizes J c for any specified value of D c 2 ; the solution is 

D = 1/2 (E wj Mj M/T y l [ Z Wi (Mj 2 - B? - < > ) Mj + D c 2 E M{\ . (6) 

We may now select the value of D that simultaneously satisfies equation (6) and the constraint 

D 2 = D c 2 . (7) 

A solution for D selected in this manner is not guaranteed to be a minimum for the loss function J. 
We do expect, however, that a vector D c chosen to satisfy equations (6) and (7) simultaneously 

will be both a good estimate for the solution that would minimize J and a good estimate for the true 
bias. In particular, we note that in the limit of noise-free data, there exists a value of D (namely, 
the true bias) that simultaneously zeroes J and J c while also satisfying the constraint of equation 

(7). Because J and J c are by definition nonnegative, a solution that zeroes them both is a minimum 
for both. It follows, therefore, that D c will be the exact solution for D in the limit of zero noise. 

With the constraint of equation (7) applied, equation (6) can be rewritten as 

D c = U + D c 2 V , (8) 

where 

U = 1/2 (E Wi Mj Mi T ) 4 E Wj (Mj 2 - B ; 2 - < ej > ) M { , (9) 

V = 1/2 (E wj Mj MT y l E Wj Mj . (10) 

Solving equation (8) for D c first requires a solution for D c 2 . This can be found as follows. 
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Multiplying each side of equation (8) with its transpose results in 


D c 2 = D c t D c = U T U + 2U t VD c 2 + V t V(D c 2 ) 2 . (11) 

This can be written as a simple quadratic in D c 2 : 


y = a (D c 2 ) 2 + b D c 2 + c = 0 , 

(12) 

where 

a = V T V , 

(13) 

b = (2 U T V - 1 ) , 

(14) 

c = U T U . 

(15) 


The solutions of equation (12) are 

D c 2 = [-b± (b 2 - 4ac) 1/2 ]/2a . (16) 

That the negatively signed solution is the proper solution can be demonstrated in the following 
manner. Any real solution for D c 2 must have b 2 > 4ac. Parameters a and c are nonnegative 

because V and U are real vectors. From this it follows that Ibl > (b 2 - 4ac ) 1/2 . In order that D c 2 
be nonnegative, it follows in turn that b < 0. Both solutions for D c 2 are therefore nonnegative. 

For cases in which neither bias nor measurement error exists, c = 0 and b = -l. The negatively 
signed solution is clearly required to give D c 2 = 0 for such cases. The parabola y(D c 2 ) defined in 
equation (12) remains the same in form (i.e., open upwards, with two nonnegative crossings of the 
D c 2 -axis) for cases of nonzero bias and error. Any physically realizable measurement situation 

may be imagined as achieved by a continuous change through a family of situations that starts with 
zero bias and zero error. During this change, the function y, and therefore the two solutions for 
equation (16), change smoothly. The negatively and positively signed solutions of equation (16) 
do not merge (and therefore cannot switch places) with each other unless the family of situations 
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contains a "point" having b 2 = 4ac. An incremental shift in the family to a neighboring one having 
b 2 > 4ac everywhere is equally acceptable and removes the conceptual possibility of sign switching 
for the true solution. The negatively signed solution is therefore always the correct solution for 
equation (16). Equations (16) and (8) can thereby be solved to obtain D c . 

The value of D c may be used either as an acceptable estimate of the magnetic bias or as an initial 

estimate in an iterative algorithm to find the value of D that minimizes the loss function J. The 
iterative algorithm is as follows. D is expressed as 

D = D 0 + d , (17) 

where D 0 is an estimate (initially set to D c ) of D, and d is a correction to this estimate. Equation 
(3) may then be written as 

ej = (Mj2 + D 0 2 - 2 MT D 0 - 2 MT d +2D Q Td+ d 2 ) - Bj 2 . (18) 

The loss function J becomes 

J=l/2Ew i [(M i 2-B i 2 +D q 2 -2MJD 0 + d 2 - <e> ) - 2 (M ; - D 0 ) T d ] 2 . (19) 

Taking the gradient of J with respect to the components of d yields 

VJ = £ Wj [ ( Mj 2 - B ; 2 + D 0 2 - 2 MT D 0 + d 2 - <e> ) - 2 (Mj - D 0 ) T d ] 

[ 2 d - 2 (Mj - D 0 ) ] . (20) 

Neglecting the terms in equation (20) that are quadratic or cubic in d yields 

VJ = 2 £ Wj { ( MT - BT + D 0 2 - 2 MT D 0 - <e ; > ) d 

-[(M.2-BT +D 0 2 -2MTD 0 -<ej>) - 2 (Mj - D 0 ) T d ] (M. - D 0 ) ) , (21) 
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which may be rearranged to read 


VJ = 2 (F d - Z ) , (22) 

where 

F = Lw i [2(M i -D 0 )(M i -D 0 ) T + (M { 2 - B ; 2 +D 0 2 -2MTD 0 -< ei >)/], (23) 

Z = E Wj [ ( Mj 2 - B; 2 + D 0 2 - 2 Mj T D 0 - <e ; > ) (M ; - D 0 ) ] , (24) 

and / is the 3-by-3 identity matrix. The desired solution for d is found by setting VJ to zero and 
multiplying both sides of the resulting equation by F' 1 : 

d = F 1 Z . (25) 

The loop is closed by adding d to D 0 , thus obtaining a new estimate for D. The process continues 
until d is made arbitrarily close to zero. As noted above, we expect that a solution for D obtained 
using this iterative technique should not differ significantly from D c . 


The uncertainty, AD, associated with the estimated value of D should be of order a/ N 1/2 , where 

a is the root-mean-square (rms) magnetometer measurement error, and N is the total number of 
measurements made. A more accurate estimate of the uncertainty of D can be made as follows. If 
equation (25) had resulted from a "standard" least-squares problem, i.e., one in which all measure- 
ment errors were associated with the vector Z while the matrix F was composed of known 

constants or model parameters, then the uncertainty AD would be given by 


AD = [2J(D)/(N-1)] 1/2 W , (26) 

with W defined as the vector with components equal to the square root of the diagonal components 
of F" 1 . If all measurement errors are "extracted" from F 1 , it reduces to a matrix X defined as 
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(27a) 


X = [Lw.2(M.-D)(M r D) 7 ]’ 1 

= [Ew i 2(A i B i )(A i B i ) T ] 1 . (27b) 

We therefore take our estimate of AD to be as given by equation (26), with W defined as the square 
root of the diagonal of X. In Section HI we give an evaluation of AD more directly comparable 
with the simple estimate a/N 1 ^ 2 . 


m. ANALYTICAL ERROR MODEL 

In this section we present a simple analytical model for the magnetometer errors and use it to 
derive expressions for the parameters w i and < e i >. We further apply this model to a simple 
observing scenario and calculate the corresponding uncertainty AD. In principle, errors associated 
with both Mj and Bj exist, the latter resulting from errors both in the geomagnetic field model and 
the spacecraft ephemeris. For simplicity, we assume here that the total error is sufficiently 
dominated by the error in M i that the error in B ; can be neglected. We define r as the error vector 

associated with and assume it to be composed of three independent random variables, each 
having Gaussian distribution with mean zero and variance a 2 . The full variance of i*j is 3a 2 . 

Using this model, the variable ej is related to r ; by the relation 

e ; = ( Mj - Dj ) 2 - Bj 2 

= (B i + r i ) 2 - Bi 2 

= 2 BT rj + rj 2 . (28) 

The expectation value of e ; is given by 
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< ej > = / G(rj) e ; dv r 

= jG(r i )(2B i Tr i + rj 2 )dv r 
= /G(r i )(2B i r i cos(0)+ r; 2 )dv r 

= 3a 2 , (29) 

where G(rj) is the Gaussian distribution function (proportional to exp(-r i 2 /2a 2 )), dv r is a volume 
element in r-space, and 0 is the angle between B ; and r ; . (The integral J f cos(0) dv r yields zero 

for any f(r) that is spherically symmetric.) The variance, s?, of ej is in turn given by 

s ; 2 = J G(r) ( e; - < e; > ) 2 dv r 

= J G(rj) [ 2 B ; rj cos(9) + r { 2 - 3c 2 ] 2 dv r 
= ! G(rj) [ 4 Bj 2 rj 2 cos 2 (0) + 2 Bj r. cos(0) (r 4 2 - 3a 2 ) 

+ r ; 4 - 2 r 2 (3a 2 ) + (3a 2 ) 2 ] dv r 
= 4/3 Bj 2 (3a 2 ) + 5 a 2 (3a 2 ) - 2 (3a 2 ) 2 + (3a 2 ) 2 

= 2 a 2 ( 2 Bj 2 + 3a 2 ) . (30) 

The weight Wj may be set equal to s ; ' 2 . We assumed for the derivation of equations (29) and (30) 

that the error in Mj dominated that in B ; . This was done for simplification and is in no way essen- 
tial. Allowing an error in Bj having components with Gaussian distribution of variance a b 2 yields 

the following expressions for < e ; > and s ; 2 : 

<ej> = 3 a 2 - 3a b 2 , (31) 

Sj 2 = 2 a 2 ( 2 B ; 2 + 3a 2 ) + 2a b 2 (2B ; 2 + 3a b 2 ) + 2(3a 2 3a fe 2 ) . (32) 
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For illustration, we present here typical values of B, a, and a b for the specific case of the HST 

mission. The rms noise level associated with the HST magnetometers has been determined to be 
approximately 0.01 gauss (Reference 2). The geomagnetic field model to be used is the Interna- 
tional Geomagnetic Reference Field (IGRF) for 1975 (Reference 3). A typical value for B at an 
altitude of 445 kilometers is about 0.3 gauss. Typical errors to be associated with the reference 
field can be estimated on the basis of the work by Trombka and Cain (Reference 4), wherein the 
IGRF for 1975 was compared with that for 1965. The rms error at 445 kilometers was found to be 
approximately 0.002 gauss. Launch date for the HST mission is currently scheduled for June 
1989, 15 years after the completion of the 1975 IGRF. Assuming linear error growth with time, 
the error to be associated with the reference field at HST launch should be about 0.003 gauss. The 
error in the estimate of the local geomagnetic field value due to spacecraft position error can be 

estimated as l(V*B)l l(AR)l, where V B is the divergence of B, and AR is the ephemeris error for 
the spacecraft. The definitive ephemeris for the HST is expected to have an rms accuracy of 0.2 
kilometer (Reference 5). Given that the near-Earth dipole field falls off approximately as R' 3 , the 

magnitude of V B can be estimated as 3 B / R 0 , where R 0 is the semimajor axis of the spacecraft 

orbit. Together these yield a field error of approximately 0.03 milligauss. The error in B ; is 
dominated by the error in the field model. The ratio o 2 / a b 2 is approximately 10, indicating that 
neglecting a fe 2 in equations (31) and (32) would be a reasonable approximation. 


We now use equations (26), (27b), and (30) to provide a simple expression for AD for an 
idealized observing scenario. We consider a situation in which the geomagnetic field is constant 
and the spacecraft changes attitude in such a way as to have the field aligned with each of the 

principal spacecraft axes for one-third of the measurements. For this scenario the variance s ; 2 is the 

same for all i, and Wj may be set to 1. The loss function J(D) and matrix X are then for large N 
given by 


J(D) = 1/2 Ns 2 = N c 2 ( 2 B 2 + 3d 2 ) , (33) 

X = (2/3NB 2 /)' 1 , (34) 
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from which it follows that each component of AD is given by 


AD = [2J(D)/N] 1/2 W 

= [ 2 a 2 ( 2 B 2 + 3a 2 ) ] 1/2 ( 2/3 N B 2 )' 1/2 

= cl N 1/2 6 1/2 [ 1 + ( 3o 2 /2 B 2 ) ] 1/2 . (35) 

The estimate specified above indicates that, for situations in which the measurement errors are small 
compared with the field strength, the algorithm described in this article will provide an estimate for 
the spacecraft bias that is about a factor of 6 1/2 (i.e., ~ 2.5) less accurate than that which could be 
obtained with the same measurements if attitude information were independently available. 


In situations where the use of the algorithm described in this article is intended as a preliminary 
calibration of the spacecraft bias, many of the details of algorithm may be neglected. As mentioned 

earlier, setting < ej > = 0 and w ; = 1 is the simplest approach in practice. If B 2 is fairly constant, 
as should be roughly true (to within a factor of about 4) for a spacecraft at constant altitude, then 
taking the weights for all variables e ; to be equal would not be a bad approximation. It would be 
inappropriate, however, for spacecraft with highly eccentric orbits. Intuitively, it might be 
expected that neglecting < ej > could introduce errors in the solution for D of perhaps as much 
as a. (Our numerical examples actually show this estimate to be too high.) Errors of this order are 
probably not significant for preliminary estimates. In the HST mission, for example, for which 
final spacecraft bias calibration will be performed using high-accuracy attitude information, the 

algorithm described here will be used with < ej > = 0 and w ; = 1 and without iteration. 


IV. NUMERICAL EXAMPLES 

In this section we present the results of a number of numerical simulations that illustrate the 
use of the algorithm and verify its accuracy. After verifying that the algorithm obtains the exact 
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solution when no noise is added to the measurements, we ran six test scenarios using noisy data. 
Each scenario comprises K bias determinations for each of two magnetometers, with each determi- 
nation based upon N data points. Both K and N have the value 100 for all of the scenarios. The 
strength of the geomagnetic field is ~ 0.35 gauss in all cases, i.e., typical of that experienced by a 
spacecraft in near-Earth orbit, and the true values for spacecraft bias at the two magnetometer 
locations are 

Dj = ( 0.005, -0.015, 0.010 ) T gauss , (36) 

D 2 = (-0.170, 0.280, 0.220 ) T gauss . (37) 

We selected these values to test the algorithm with bias values substantially below and approxi- 
mately equal to the ambient field strength. For each scenario, we calculated the average error in 
the estimated bias (for both D c and D) and the standard deviation of these errors for each magne- 
tometer. For calculating D, we defined convergence as Idl < 0.001 <r/N 1/2 . 

We divide the six scenarios into two groups. The scenarios in the first group (scenarios 1, 2, 
and 3) all use a constant geomagnetic field. We are therefore able to compare the results of these 
scenarios directly with the analytic calculations presented in Section III. The scenarios in the 
second group (scenarios 4, 5, and 6) use a geomagnetic field constructed to simulate that which 
could be experienced by a spacecraft in near-Earth orbit. The height of the simulated orbit is 
522 kilometers; the inclination is 28.5 degrees; and the right ascension of ascending node, 
argument of perigee, and latitude at epoch time (i.e., time of the first data point) are all 0 degrees. 
The geomagnetic field experienced by a spacecraft with this orbit is reasonably simulated (at least 
for our purposes) using the following parameterization: 


j j = 0.01+ 0.17 cos(0.) gauss , 

(38) 

2 i = -0.19+ 0.15 sin(0.) gauss, 

(39) 

>3 j = 0.20 + 0.07 sin(0.) gauss , 

(40) 


where 0. = 7.2 ( i - 1) degrees. The sample spacing corresponds roughly to one data point per 
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minute over a single orbit We present details for the individual scenarios in the following para- 
graphs and the principal numerical results in Table 1 (for scenarios 1, 2, and 3) and Table 2 (for 
scenarios 4, 5, and 6) at the end of this section. The specific entries in the tables are average error, 
standard deviation of the errors, estimated uncertainty based on equation (26), and for scenarios 1 
and 2 the "theoretical" uncertainty calculated using equation (35). 

Scenario 1 . This scenario is a numerical execution of the idealized case described in Section HI. 
It is designed to demonstrate the algorithm's accuracy for situations where the magnetometer 

measurements are uniformly distributed. The values for the parameters B and a are 0.35 gauss 
and 0.01 gauss, respectively. As can be seen from Table 1, the magnetic bias estimates are equally 
accurate for both magnetometers, and no significant difference exists between the estimates D and 
D c . Typically, only two iterations are required to obtain the estimate D after D c has been deter- 
mined. The standard deviation for each of the error components is ~ 0.0018 gauss, in contrast to 
the "theoretical" value of ~ 0.0025 gauss predicted by both equations (26) and (35). This differ- 
ence may be attributable to the qualitative nature of the concept of "extracting" the measurement 
errors from the matrix F' 1 that was used in deriving equation (27b). The average error in the 
estimated bias for each component of D (or D c ) is zero to within plus or minus about twice 
1.8 a / (N K) 1/2 (i.e., ~ twice 0.00018 gauss), which is what would be expected given the 
standard deviation of the errors. 

Scenario 2 . This scenario is a repetition of scenario 1, except that in this case the parameter 
< e } > is set to zero. This scenario is designed to determine whether the algorithm with this 

simplification provides results of essentially the same accuracy as does the full algorithm. This, in 
fact, proves to be the case. The standard deviation of the errors is essentially the same as for 

scenario 1 (i.e., ~ 1.8 a / N 1 ^ 2 ), and the average of the errors for any given component of D or D c 

is shifted only slightly (~ 0.0006 gauss) from zero. This offset, though statistically significant, is 
small compared with the typical uncertainty of ~ 0.002 gauss for any single determination. It 

appears that a user of the algorithm can apply the simplification of neglecting both < e ; > and the 

steps of iterating between D c and D without substantially corrupting the results of the calculation. 



Scenario 3 . This scenario is designed to investigate the observability problems that the algorithm 
encounters if an observing run with insufficient variation of the geomagnetic field is used. The 
field strength and noise levels are the same as for scenario 1, but the three attitudes used are all 
within 14 degrees of one another. Specifically, the right ascension, declination, and roll for the 
three attitudes are, in degrees, (0, 0, 0), (10, 0, 0), and (0, 10, 0); and the geomagnetic field is 
directed toward sky coordinates (0, 0). As shown in Table 1, component 1 of the standard 
deviation vector for the errors of D is substantially smaller (by a factor of between 7 and 9) than 
components 2 and 3. Furthermore, the same kind of ratio between components applies for the 
uncertainty estimate derived using equation (26). Despite this basic agreement (i.e., propor- 
tionality) between the standard deviation vector and the uncertainty estimate, the average errors of 
the components of D are larger in size (i.e., offset from zero) by amounts that are statistically 
significant when compared with their expected rms size of K' 1/2 times the standard deviation. 

This is true in particular for component 1 for magnetometer 1 and components 2 and 3 for 
magnetometer 2. Although the offsets are statistically significant, they are all still smaller than the 
corresponding standard for the given component of D. The calculation of D remains basically 
good, but significantly degraded, for this scenario. In contrast to these results for D, the bias 
estimate D c is a poor estimate of the true bias (particularly in the case of magnetometer 2, for which 
the true bias is large). For the scenario presented here, the averages of the errors for components 
2 and 3 of D c for magnetometer 2 are nearly 10 times as large as the rms error of a single measure- 
ment. These results emphasize the need for the user to arrange for data acquisition in a manner that 
guarantees substantial field variation in the spacecraft’s frame. They further indicate that if such 
uniform data acquisition is for some reason not possible, then application of the full algorithm (i.e., 
continuing on to calculate D after D c ) is imperative. 

Scenario 4 . This scenario, as well as the next two, is designed to test the algorithm using the 
more realistic geomagnetic field model described at the beginning of this section (i.e., as specified 

by equations (38) through (40) ). The parameters Wj and < e ; > are set in accordance with equa- 
tions (29) and (30), with a again equal to 0.01 gauss. As can be seen in Table 2, the results for 
scenario 4 with respect to accuracy of the solution are essentially the same as those for scenario 1 
(component 3 of the error for both D and D c for magnetometer 2 is just slightly greater than three 
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times the expected rms value). These results demonstrate that the algorithm works well for realistic 
ambient magnetic field values, as well as for the idealized cases such as scenario 1. 

Scenario 5 . This scenario is a repetition of scenario 4, but with w ; set to 1 and < ej > set to 0. 

This scenario is designed with the same purpose as scenario 2, i.e., to check the accuracy of the 
algorithm when it is applied with the indicated simplifications. Comparing the results of this 
scenario and those of scenario 4 indicates (as was noted for scenario 2) that a mild degradation of 
results occurs when the simplified algorithm is used, but that this degradation is quite small. In 
actual operations, use of the simplified algorithm would almost certainly be acceptable. For the 

realistic error estimate being used in these scenarios (i.e., a ~ 0.01 gauss), the final uncertainty in 
the bias estimate (using either D or D c ) is of order 0.003 gauss. 

Scenario 6 . This scenario is a repetition of scenario 4, but with a set to 0.1 gauss. The purpose 
is to test the robustness of the algorithm when used with very noisy data; the noise level used is 
about one-third of the ambient field strength. The algorithm requires approximately six iterations 
to converge and yields estimates for the bias with errors of order 0.04 gauss. In one case the 
algorithm failed to converge even after 10 iterations. As is clear in Table 2, the components of the 
vector of average errors are offset from zero by as much as 10 times the expected rms value. There 
is no indication that the estimate D is any better determined than is D c . 

The basic conclusions to be the drawn from the six scenarios are that (1) the algorithm works 
very well for cases in which magnetometer noise levels are like those expected for actual spacecraft 

operations, (2) results obtained using the simplifications of setting w } to 1 and < e ; > to 0 and not 

iterating beyond the first estimate D c are very nearly as good as those obtained using the full 

algorithm, and (3) the user of the algorithm should make every effort to arrange that magnetometer 
data acquisition be done in a way that provides substantial variation in all components. The results 
of scenarios 4 and 5 indicate that the kind of magnetic variation that a spacecraft with constant 
attitude in near-Earth orbit automatically experiences over a single orbit is adequate to the needs of 
the algorithm. 
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Table 1. Statistical Results from Scenarios 1 through 3 


Magnetometer: 


1 



2 


Component: 

1 

2 

3 

1 

2 

3 

Scenario 1 







Average Error 
D: Standard Deviation 

Estimated Uncertainty 

0.00038 

0.00179 

0.00242 

0.00038 

0.00184 

0.00242 

0.00000 

0.00182 

0.00239 

0.00008 

0.00166 

0.00246 

0.00033 

0.00190 

0.00246 

0.00027 

0.00181 

0.00242 

Average Error 
D c : Standard Deviation 

0.00037 

0.00178 

0.00042 

0.00183 

-0.00002 

0.00181 

0.00029 

0.00165 

-0.00001 

0.00188 

-0.00001 

0.00182 

Estimated Uncertainty 

0.00242 

0.00242 

0.00239 

0.00246 

0.00246 

0.00242 

Theoretical Uncertainty 

0.00246 

0.00246 

0.00243 

0.00246 

0.00246 

0.00243 

Scenario 2 







Average Error 
D: Standard Deviation 

Estimated Uncertainty 

0.00081 

0.00179 

0.00242 

0.00081 

0.00184 

0.00242 

0.00043 

0.00182 

0.00238 

0.00051 

0.00166 

0.00245 

0.00076 

0.00190 

0.00245 

0.00070 

0.00181 

0.00242 

Average Error 
D c : Standard Deviation 

0.00080 

0.00178 

0.00085 

0.00183 

0.00041 

0.00181 

0.00072 

0.00165 

0.00041 

0.00188 

0.00042 

0.00181 

Estimated Uncertainty 

0.00242 

0.00242 

0.00238 

0.00245 

0.00245 

0.00242 

Theoretical Uncertainty 

0.00246 

0.00246 

0.00243 

0.00246 

0.00246 

0.00243 

Scenario 3 







Average Error 
D: Standard Deviation 

Estimated Uncertainty 

0.00093 

0.00209 

0.00242 

0.00000 

0.01530 

0.01958 

-0.00127 

0.01677 

0.01943 

0.00018 

0.00149 

0.00235 

-0.00807 

0.01342 

0.01899 

-0.00955 

0.01386 

0.01885 

Average Error 
D c : Standard Deviation 

0.00077 

0.00178 

0.00146 

0.01265 

-0.00150 

0.01405 

0.01664 

0.00552 

-0.09821 

0.01369 

-0.09171 

0.01398 

Estimated Uncertainty 

0.00242 

0.01959 

0.01944 

0.00295 

0.02386 

0.02368 


Note: All values are in gauss. 
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Table 2. Statistical Results from Scenarios 4 through 6 


Magnetometer: 


1 



2 


Component: 

1 

2 

3 

1 

2 

3 

Scenario 4 







D: 

Average Error 
Standard Deviation 
Estimated Uncertainty 

-0.00003 

0.00272 

0.00375 

0.00035 

0.00240 

0.00317 

0.00032 

0.00212 

0.00292 

0.00047 

0.00285 

0.00373 

-0.00023 

0.00243 

0.00315 

0.00077 

0.00220 

0.00290 


Average Error 
Standard Deviation 

-0.00006 

0.00270 

0.00045 

0.00239 

0.00023 

0.00211 

0.00097 

0.00282 

-0.00050 

0.00241 

0.00076 

0.00219 

Estimated Uncertainty 

0.00375 

0.00317 

0.00292 

0.00373 

0.00315 

0.00290 

Scenario 5 







D: 

Average Error 
Standard Deviation 
Estimated Uncertainty 

-0.00008 

0.00270 

0.00378 

0.00056 

0.00243 

0.00306 

0.00087 

0.00216 

0.00323 

0.00054 

0.00285 

0.00376 

0.00003 

0.00244 

0.00304 

0.00128 

0.00220 

0.00321 

»c- 

Average Error 
Standard Deviation 

-0.00011 

0.00269 

0.00066 

0.00241 

0.00077 

0.00215 

0.00102 

0.00283 

-0.00016 

0.00242 

0.00122 

0.00219 

Estimated Uncertainty 

0.00378 

0.00306 

0.00323 

0.00376 

0.00304 

0.00321 

Scenario 6 







D: 

Average Error 
Standard Deviation 
Estimated Uncertainty 

0.00086 

0.03783 

0.03798 

0.02297 

0.03083 

0.03188 

0.04112 

0.03081 

0.02992 

0.00856 

0.04478 

0.03749 

0.01891 

0.03685 

0.03146 

0.04873 

0.04229 

0.02953 


Average Error 
Standard Deviation 

-0.00252 

0.02238 

0.02606 

0.01896 

0.02527 

0.01700 

0.03283 

0.02563 

-0.00035 

0.02270 

0.02891 

0.02359 

Estimated Uncertainty 

0.03817 

0.03203 

0.03007 

0.03847 

0.03228 

0.03031 


Note: All values are in gauss. 
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V. CONCLUSIONS 


We have presented a simple algorithm for the in-flight determination of the magnetic bias of a 
spacecraft. The algorithm, developed for use during the HST mission, is independent of any 
attitude estimates and requires no telemetry from spacecraft sensors other than the magnetometers 
at the locations where the bias is being determined. We have presented an analytical model for the 
errors in a typical observing scenario, as well as the results from a number of numerical studies 
demonstrating the general stability and accuracy of the algorithm. The algorithm works well and 
should be a useful tool for preliminary spacecraft magnetic bias determination for a variety of 
spacecraft missions. 

The work reported in this article was supported in part by NASA contract NAS 5-26685 for the 
development of ground support software for the Hubble Space Telescope mission. 
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